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ABSTRACT 

Using the RXTE/PCA data we study the fast variabiHty of the reflected emis- 
sion in the soft spectral state of Cyg X-1 by means of Fourier frequency resolved 
spectroscopy. We find that the rms amplitude of variations of the reflected emission 
has the same frequency dependence as the primary radiation down to time scales of 
^ 30 — 50 msec. This might indicate that the reflected flux reproduces, with nearly flat 
response, variations of the primary emission. Such behavior differs notably from the 
hard spectral state, in which variations of the reflected flux are significantly suppressed 
in comparison with the primary emission, on time scales shorter than ~ 0.5 — 1 sec. 

If related to the finite light crossing time of the reflector, these results suggest that 
the characteristic size of the reflector - presumably an optically thick accretion disk, 
in the hard spectral state is larger by a factor of <; 5 — 10 than in the soft spectral 
state. Modeling the transfer function of the disk, we estimate the inner radius of the 
accretion disk ~ lOORg in the hard and Rin ^ lORg in the soft state for a IOMq 
black hole. 

Key v^rords: accretion, accretion disks - black hole physics - stars: binaries: general 
- stars: individual (Cygnus X-1) - X-rays: general - X-rays: stars 



1 INTRODUCTION 

, The importance of the reprocessed/reflected component in 
' the X-ray spectra of accreting X-ray sources for exploring 
the geometry of the accretion flow is well known. Reflec- 
tion of the primary Comptonized radiation from neutral or 
partially ionized matter located in the vicinity of the com- 
pact object - presumably the optically thick accretion disk, 
leads to appearance of characteristic features in the spec- 
tra of X-ray binaries. The main signatures of the emission 
reflected from cold neutral medium are well known - the 
fluorescent Kc line of iron at 6.4 keV, iron K-edge at 7.1 
keV and a broad hump at ~ 20 — 30 keV (Basko, Sunyaev 
& Titarchuk 1974, George & Fabian 1991). The exact shape 
of these spectral features in the X-ray binaries depends on 
ionization state of the reflecting medium and might be mod- 
ified by strong gravity effects and intrinsic motions in the 
reflector (e.g. Fabian et al., 1989). The amplitude of the re- 
flection signatures depends primarily on the ionization state 
and the solid angle subtended by the reflector as seen from 
the source of the primary radiation. 

Recently, Revnivtsev, Gilfanov & Churazov (1999, here- 
after Paper I) proposed Fourier frequency resolved spectral 



analysis to study spectral variability in the X-ray binaries. 
Although interpretation of the Fourier frequency resolved 
spectra in general is not straightforward and requires a pri- 
ori assumptions to be made, one of the areas where this 
method can be efficiently used is the fast variability of the 
reflected component. It has been found that in the low spec- 
tral state of Cyg X-1 and GX339-4 the energy spectra cor- 
responding to the shorter time scales {^0.1 — 1 sec) show 
less reflection than those of longer time scales (Revnivtsev, 
Gilfanov & Churazov 1999, 2000). The simplest, although 
not unique, interpretation of this result is smearing of the 
short term variations of the reflected emission due to finite 
light crossing time of the reflector. Based on the Fourier fre- 
quency dependence of the equivalent width of the Fe Ka 
fluorescent line Revnivtsev et al. (1999) estimated the char- 
acteristic size of the reflector: ~ 80 — IGORg for a 20 — IOMq 
black hole. 

In this paper we investigate the fast variability of the 
reflected component in the soft spectral state of Cyg X-1 
and compare it with that in the hard spectral state. 
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Figure 1. The spectral energy distribution of Cyg X-1 in the 
soft (filled circles) and hard (open circles) spectral state. The data 
of nearly simultaneous ASCA and RXTE observations on March 
26, 1996 (hard state) and May 30. 1996 Isoft state) . A source 



distance of 2 kpc was assumed (Gierlinski et al. 1999) 



2 OBSERVATIONS AND DATA ANALYSIS. 

We used publicly available data of Cyg X-1 observations 
with the Proportiona l Counter Array aboard the Ros si X- 
ray Timing Explorer ( Brandt Rotschild fc Swank 1996[ ) per- 



formed in June, 1996 during the soft spectral state of the 
source. The list of the observations is given in Table |l[ The 
total live time was ~ 11.5 ksec. The "Generic Binned" mode 
data in configuration B_4ms_8A_0_35_H, with time resolu- 
tion of ~ 4 msec and covering 2.9-13.1 keV energy range 
was used for frequency resolved spectral analysis. 

The data screening and selection was performed using 
FTOOLS 4.2 with standard screening criteria recommended 
by RXTE GOF. The frequency resolved spectra were cal- 
culated following the prescription detailed in Paper I. The 
dead time corrected value of the white noise level was de- 
termined from fitting of the power spectra in the 300-1000 
Hz frequency ra nge. The spect ral analysis was performed in 
XSPEC v.10.0 ( lArnaud 199^ ) with version 3.5 of the PCA 
response matrix. A uniform systematic error of 0.5% was 
added quadratically to the statistical error in each energy 
channel. 

The observations of the source during the hard state 
discussed in the text were performed between March 26 and 
31, 1996. The details of these observations are given in Paper 
1. 



3 RESULTS. 

The broad band energy spectra of Cyg X-1 in the hard and 
soft spectral state are shown in Fig.bJ. The spectra were ob- 



Figure 2. The power density spectra of Cyg X-1 in the soft 
(filled circles) and hard (open circles) spectral state in the 3- 
13 keV energy band. The power spectra are plotted as (power 
density) X frequency, i.e. in units of rms^/Hz X Hz and are cut 
at approximately half of the Nyquist frequency (32 and 128 Hz for 
the hard and soft state respectively). No correction for binning 
and aliasing effects has been applied. 



Table 1. The list of RXTE observations used for the analysis. 



Obs.lD 


Date (UT) 


Start 


End 


Expos., s* 


10512-01-05-00 


04/06/96 


20:49:36 


21:43:44 


1249 


10512-01-07-00 


16/06/96 


00:07:44 


00:21:04 


643 


10512-01-07-02 


16/06/96 


04:55:44 


05:17:04 


1113 


10512-01-08-01 


17/06/96 


01:44:32 


01:56:00 


532 


10512-01-08-02 


17/06/96 


04:56:32 


05:17:52 


1081 


10512-01-08-00 


17/06/96 


08:08:32 


08:44:00 


1871 


10512-01-09-02 


18/06/96 


03:21:36 


03:36:00 


675 


10512-01-09-00 


18/06/96 


06:34:40 


07:01:52 


1461 


10512-01-09-01 


18/06/96 


09:46:40 


10:26:56 


2170 



* Dead time corrected PCA exposure time. 

tained using the data of overlapping ASCA and RXTE ob- 
servations of the source on March 26, 1996 (hard state) and 
on May 30, 1996 (beginning of the soft state). 

The power spectra of Cyg X-1 in the 3-13 keV energy 
band in the hard and soft spectral states obtained from 
the complete sets of the data used for the frequency re- 
solved spectral analysis in this paper and in Paper I are 
shown in FigJ^. The power density is plotted in units of 
power X frequency representing squared fractional rms at a 
given frequency per factor ~ e in frequency. This way of 
representing the power spectra most clearly characterizes 
relative contribution of variations at different frequencies to 
the total observed rms. 

The Fourier frequency resolved spectra for the soft spec- 
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Table 2. Dependence of the best fit parameters upon Fourier frequency 



Frequency 




soft state 






hard state 




range, Hz 


phot.ind. 


EW2, eV 


EW/EW(< lHz)3 


phot.ind. F^ 


EW^, eV 


EW/EW(< lHz)3 


0.016-0.078 


2.60 ±0.032 


494 ± 83 


1.05 ±0.18 


1.97 ±0.006 


145 ± 16.6 


0.97 ±0.11 


078-0 20 


2.60 ± 0.022 


484 ± 61 


1 03 + 1 3 


1 95 + 004 


153 + 130 


1 02 + 09 


0.20-0.45 


2.63 ±0.018 


482 ± 48 


1.02 ±0.10 


1.94 ±0.004 


152 ± 11.3 


1.02 + 0.08 


0.45-1.0 


2.62 ±0.014 


424 ± 38 


0.90 ± 0.08 


1.94 ±0.004 


149 ± 10.4 


1.00 ±0.07 


1-2 


2.63 ±0.012 


445 ± 33 


0.94 ± 0.07 


1.89 ± 0.004 


121 ± 9.8 


0.81 ±0.07 


2-4 


2.64 ±0.011 


429 ± 31 


0.91 ±0.07 


1.86 ± 0.004 


121 ± 9.9 


0.81 ±0.07 


4-8 


2.64 ±0.011 


413 ± 31 


0.88 ±0.07 


1.81 ± 0.004 


107 ±11.4 


0.71 ±0.08 


8-16 


2.58 ±0.011 


462 ± 35 


0.98 ±0.07 


1.70 ±0.007 


75 ±16.6 


0.50 ±0.11 


16-32 


2.52 ±0.018 


410 ± 59 


0.87 ±0.13 


1.69 ± 0.017 


47 ±37.8 


0.31 ±0.25 


32-128 


2.43 ± 0.050 


233 ± 172 


0.49 ±0.37 









The errors are la for one parameter of interest. 

^ — the power law photon index; 

^ - the equivalent width of the 6.4 keV line, eV; 

^ - EW(< IHz) - equivalent width averaged in the 0.016-1.0 Hz frequency range 



tral state were obtained in 10 frequency bins of logarithmi- 
cally equal width in the 0.016-128 Hz frequency range. In 
order to study the frequency dependence of the amplitude 
of the reflected component we fit the spectra in the 3-13 
keV band with a simplified model consisting of an absorbed 
power laviQwith superimposed Gaussian line at 6.4 keV. The 
low energy absorption was fixed at Nh = 6 • 10^^ cm~^. The 
line width (standard deviation for a Gaussian profile) was 
fixed at 1 keV which corresponds to the average value for 
the high state observations used for the analysis. The power 
law photon index and line fiux were the only free parameters 
of the fit. To facilitate comparison with the hard state data 
we reanalyzed the data set used in Paper I in identical fre- 
quency bins covering 0.016-32 Hz frequency range (the time 
resolution of the hard state data was ^ 16 msec) and using 
the same spectral model. The width of the line for the hard 
state spectr a was fixed at the corresponding averag e value 
of 0.8 keV (iGilfanov, Revnivtsev fc Churazov 1999|). 



Such a spectral model is obviously oversimplified. Nei- 
ther it is justified from the physical point of view. It is, 
however, suitable to quantify the amplitude of the charac- 
teristic "wiggle" usually seen in the spectra of the accreting 
X-ray sources in the ~ 5 — 15 keV energy range and com- 
monly attributed to the effects of refiection. Use of more 
sophisticated models is restricted by insufficient statistics 
(especially in the soft spectral state) and the low number 
of energy channels below ~ 13 keV in the B_4ms_8A_0_35_H 
configuration used in the most of the soft state observations. 

The best fit values of the spectral parameters for both 
spectral states are given in Table ^ The dependence of the 
equivalent width on the Fourier frequency is shown in Fig.^ 
Variations of the parameters of the spectral model, in par- 
ticular the change of the line centroid from 6.0 to 6.7 keV 
and reducing the intrinsic line width to zero, do not change 
the general trend. These variations, however, affect the par- 



* Contrary to the average spectra of Cyg X-1 in the soft state 
(Fig.hh , the contribution of the soft component to the frequency 
resolved spectra is negligible (to be discussed in more detail in a 
separate paper) . Therefore use of a power law to model continuum 
emission in the soft state is justified. 
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Figure 3. Equivalent width of the Fe fluorescent line vs. Fourier 
frequency in the soft (filled circles) and hard (open circles) spec- 
tral state. The equivalent width is normalized to the 0.01-1 Hz 
average. 



ticular values of the equivalent width and, to lesser extent, 
the shape of the curve in Fig. |^. 

According to the x'^~test two distributions, shown in 
Fig.| differ at the confidence of level of ^ 98.7% (x^ = 19.4 
for 8 d.o.f. in the 0.016-32 Hz frequency range). It should be 
noted, however, that the errorbars assigned to the Fourier 
frequency resolved spectra were propagated from the corre- 
sponding power density spectra and are likely to be some- 
what overestimated, especially in the low frequency bins 
(this fact can be noticed in Fig.^). The confidence level, 
calculated using 0.45-32 Hz frequency range is ~ 99.7% 
(X^ = 18.1 for 5 d.o.f.). 
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4 FREQUENCY DEPENDENCE OF THE 
EQUIVALENT WIDTH AND TIME 
RESPONSE OF THE REFLECTOR. 

The geometry, element abundances and ionization state of 
the reflector being fixed, the equivalent width of the Fe flu- 
orescent line determined from a conventional energy spec- 
trum is proportional to the relative amplitude of the re- 
flected component and approximately measures the solid 
angle subtended by the reflector. The equivalent width of 
the fluorescent line determined from a Fourier frequency re- 
solved spectrum measures the ratio of the rms amplitudes of 
variations of the reflected component and the primary emis- 
sion in a given Fourier frequency range. The constancy of the 
equivalent width of the line at Fourier frequencies f ^ 30 Hz 
observed in the soft state (Fig.^ implies that the rms ampli- 
tude of variations of the reflected component has the same 
frequency dependence as that of the primary radiation. This 
would naturally appear if the reflected emission was repro- 
ducing variations of the primary radiation down to the time 
scales of ~ 30 — 50 msec. Such behavior is in contrast to 
the hard spectral state, in which variations of the reflected 
flux are notably suppressed in comparison with the primary 
emission on the time scales shorter than ~ 500 msec. 

The most straightforward explanation of this effect 
would be in terms of a finite light crossing time of the reflec- 
tor Trcfl ~ /rcfl/c duc to its finite spatial extend Zrcfl. In this 
case the frequency dependence of the fluorescent line equiv- 
alent width, EW{f), is determined by the geometry of the 
primary source and the reflector. The characteristic width of 
the EW{f) is mainly defined by the spatial extent of the part 
of the accretion disk giving the main contribution to the re- 
flected emission. Based on these arguments Revnivtsev, Gil- 
fanov & Churazov (1999) estimated the characteristic size 
of the reflector in the low spectral state as Irea ~ 5 x 10* cm 
which would correspond to ~ 150-Rg for a IOMq black hole. 

Below we consider this problem in a more quantitative 
way. The time dependence of the reflected emission is deflned 
by the following relation: 

oo 

Frefl(t) = J Foit-T)T(r)dT 


where Fo{t) and Frefi{t) are primary and reflected flux, T{t) 
- the transfer function of the reflector, defined by the geom- 
etry. This function accounts for the propagation time of the 
photons from the primary source to different parts of the 
reflector and then to the observer. The Fourier transform of 
the reflected flux is: 

Fr.fiif) ^ Foif) X fif) 

where Fq, Frefi and T are corresponding Fourier transforms 
and / is Fourier frequency. The equivalent width of the fluo- 
rescent line determined from the Fourier frequency resolved 
spectra is 



EW{f) cx 



\Frefl{f)\ 

\Foif)\ 



\nf)\ 



h above a flat disk with the inner radius Rin and incli- 
nation angle i for different values of Rin and i. A Lam- 
bert law for the angular dependence of the reffected flux 
has been assumed. No general or special relativity effects 
have been taken into account. The characteristic width of 
the EW{f) dependence is mainly deflned by the distance 
+ ^in between the primary source and the inner 
edge of the disk and depends only weakly on the inclina- 
tion angle i. A small offset of the primary source, A <^ d 
or non-zero opening angle of the disk do not signiflcantly 
affect the characteristic width of Fourier transform of the 
response function (Fig.^. However, the transfer function it- 
self and the high frequency part of its Fourier transform are 
sensitive to the details of the geometry (Fig. |^ and ^ . 

In Fig.^ we compare Fourier transform of the transfer 
function of a fiat disk with inclination i — 50°, appropriate 



for Cyg X-1 (Gies & Bolton 1986), with the frequency de- 



i.e. is proportional to the Fourier transform of the transfer 
function of the reflector. 

Fig.^ shows the transfer function and its Fourier trans- 
form for an isotropic point source located at the height 



pendence of the equivalent width observed in the soft and 
hard spectral states of Cyg X-1. As seen from Fig.^ sup- 
pression of the high frequency variations in the reflected 
emission observed in the hard state can be satisfactorily 
described by reflection from a disk with an inner radius 
of Rin ~ 100-Rg around a IOMq black hole. Signiflcantly 
larger, Rin ^ 1000i?g, or smaller values of the inner radius, 
Rin ~ WRg, are inconsistent with the data. The soft state 
data, on the other hand, requires much smaller values of the 
inner radius of the disk, Rin ~ 10-Rg. 

ft should be noted that since we use the equivalent 
width of the Fe fluorescent line as a measure of the ampli- 
tude of the reflected emission, the results might be somewhat 
affected by the non-uniformity of the ionization state of the 
accretion disk with radius and geometrical effects (e.g. radial 
and azimuthal dependence of the reflection angle). Results 
of more detailed modeling of the disk transfer function and 
rigorous comparison with the data will be published else- 
where (a paper in preparation). 



5 DISCUSSION. 

Study of fast variability of the reflected emission by means 
of Fourier frequency resolved spectroscopy offers a new in- 
dependent method to probe the geometry of the accretion 
flow which complements conventional X-ray spectral analy- 
sis. The presence of a luminous soft component dominating 
the X-ray spectrum in the soft spectral state of black hole 
candidates (Fig.|l|) strongly favors small values of the inner 
radius of the disk (e.g. Gierlinski et al., 1999), in good agree- 
ment with the above result, Rin ^ lORg. In the hard spec- 
tral state conventional estimates of the inner radius of the 
accretion disk range from several tens Rg in coronal model 
(e.g. Poutanen, Krolik & Ryde, 1997; Done & Zicky, 1999) to 
several hundred Rg in ADAF model (e.g. Esin et al., 1998). 
Our result, Rin 100-Rg, falls in the middle of this range. 

The inner radius of the accretion disk determined in the 
above analysis refers to the inner radius of the "reflective" 
part of the disk, where the ionization state is such that the 
disk is capable to produce a fluorescent iron line. Therefore 
substantial change of the ionization state of the surface layer 
of the inner disk (e.g. Young et al., 1999) may have similar 
effect on the frequency dependence of the equivalent width 
of the iron line as physical change of the inner disk radius. 
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Figure 4. Transfer function (left) of the disk and its Fourier transform (right) for an isotropic point source of primary radiation located 
at h = lORg above a flat disk with inner radius of Ri„ = 10, 100 and 1000-Rg and inclination angle of 0°, 45° and 80° (not shown in the 
left panel). No relativistic effects have been taken into account. The linear quantities are given in the units of gravitational radii for a 
IOMq black hole. 




Figure 5. Dependence of the transfer function and its Fourier transform upon offset of the source of the primary radiation from the disk 
axis (offset=10Rg , flat disk - dotted line) and opening angle of the disk (offset=0, H/ R = 0.1 - dashed line and H/R = 0.5 - dash-dotted 
line) . Other parameters of the model are: inner disk radius Rin = lOORg , M = IOMq , inclination 45° . No relativistic effects have been 
taken into account. 
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Figure 6. Equivalent width of the Fc fluorescent line vs. Fourier 
frequency. Comparison of the data with the model. The data 
points are the same as in Fig.^ The model curves were calcu- 
lated for an isotropic point source at height h = lORg on the axis 
of a flat disk with inner radius of 10, 100 and 1000-Rg (assuming 
a IOMq black hole) and with inclination angle of 50° . 

In the simplified geometry of an isotropic point source 
above flat disk solid angle subtended by the reflector is ^ = 
— , ^ (in the notation of the previous section). For the 

parameters from Fig.^ respective values of the solid angle 
are: ^ 0.1 for the hard state (h = WRg, Ri„ = WORg) 
and ^ « 0.7 for the soft state (h = lORg, Rin = lOi^g). 
Equivalent width of the iron fluorescent line expected for 
these values of the solid angle and solar abundance of iron 
are several times smaller than those given in the Table ^ 
especially for the hard spectral state. However, as was noted 
above, the absolute values of the equivalent width quoted in 
Table |^ are subject to some uncertainty due to simplified 
model used for the spectral fits. Comparison of the model 
predictions with more accurately determined values of the 
equivalent width and amplitude of the refiected component 
can further constrain geometry of the accretion flow. 

Finally, we should note that interpretation of the fre- 
quency dependence of the equivalent width of the fluorescent 
line in terms of the finite light crossing time of the reflector 
is not unique. An alternative explanation might be that the 
short time scale, ^ 50-100 msec, variations appear in a ge- 
ometrically different, likely inner, part of the accretion flow 
and give a rise to significantly weaker, if any, reflected emis- 
sion than the longer time scale events presumably originat- 
ing in the outer regions. This might be caused, for instance, 
by a smaller solid angle of the refiector as seen by the short 
time scale events and/or due to screening of the reflector 
from the short time scale events by the outer parts of the 
accretion flow. However, independent of the nature of the 
fall off of the equivalent width at high frequency, the flat re- 



sponse of the reflected emission to variations of the primary 
radiation observed at low frequencies puts an upper limit on 
the spatial extent of the reflector, i.e. on the inner radius of 
the accretion disk. In particular, large values of Rin, signif- 
icantly exceeding ~ lOORg in the hard and ~ lOiig in the 
soft spectral state are excluded by our analysis. 



6 CONCLUSIONS. 

We have exploited Fourier frequency-resolved spectral anal- 
ysis to study fast variability of the reflected emission on time 
scales of ~ 100 sec - 10 msec in the soft and hard spectral 
states of Cyg X-1. Our conclusions are; 

(i) In the soft spectral state variations of the reflected 
component have the same frequency dependence of the rms 
amplitude as the primary emission up to the frequencies ^ 
30 Hz. This would be expected if, for instance, the reflected 
flux was reproducing, with flat response, variations of the 
primary radiation down to the time scales of ~ 30 — 50 
msec. The sensitivity of the present analysis is insufficient 
to study shorter time scales. 

(ii) In the hard spectral state variability of the reflected 
flux is significantly suppressed in comparison with the direct 
emission on the time scales shorter than ~ 0.5 — 1 sec (see 
also Paper I). 

(iii) Assuming that suppression of the short-term vari- 
ability of the reflected emission is caused by the finite light- 
crossing time of the reflector, we estimated the inner radius 
of the accretion disk Ri„ ~ lOOiig in the hard spectral state 
and Rin lORg in the soft spectral state. 
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